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Abstract 



We investigate models of future finite-time singularities in f(T) theory, where T is the torsion 
^^. scalar. The algebraic function f(T) is put as the teleparallel term T plus an arbitrary function g{T). 

< I ■ A suitable expression of the Hubble parameter is assumed and constraints are imposed in order to 

provide an expanding universe. Two parameters /3 and Hs that appear in the Hubble parameter are 



relevant in specifying the types of singularities. Differential equations of g{T) are established and 



^ ' solved, leading to the algebraic f{T) models for each type of future finite time singularity. Moreover, 

in ■ 

we take into account the viscosity in the fluid and discuss three interesting cases: constant viscosity, 
viscosity proportional to \/ —T and the general one where the viscosity is proportional to (—T)"", 



where n is a natural number. We see that for the first and second cases, in general, the singularities 
^— ^ ' are robust against the viscous fluid, while for the general case, the Big Rip and the Big Freeze can 

CN ■ 

be avoided from the effects of the viscosity for some values of n. 



Pacs numbers: 

1 Introduction 

Recent observations of type la supernova (SNIa) and WMAP [11 [2] indicate that our universe is 
currently undergoing an accelerating expansion, which confront the fundamental theories with great 
challenges and also make the researches on this problem a major endeavour in modern astrophysics and 
cosmology. Missing energy density - with negative pressure - responsible for this expansion has been 
dubbed dark energy. Wide range of scenarios have been proposed to explain this acceleration while most 
of them can not explain all the features of universe or they have so many parameters that makes them 
difficult to fit. In this direction we can consider theories of modified gravity i3j, or field models of dark 
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energy. The field models that have been discussed widely in the literature consider a cosmological constant 
[4] , a canonical scalar field (quintessence) J.5 , a phantom field, that is a scalar field with a negative sign of 
the kinetic term [6l [7] , or the combination of quintessence and phantom in a unified model named quintom 
[8] . In the other hand modified models of gravity provides the natural gravitational alternative for dark 
energy [9]. Moreover, modified gravity present natural unification of the early time infiation and late-time 
acceleration thanks to different role of gravitational terms relevant at small and at large curvature. Also 
modified gravity may naturally describe the transition from non-phantom phase to phantom one without 
necessity to introduce the exotic matter. Among these theories, scalar-tensor theories [10], f{R) gravity 
[11) . DGP braneworld gravity [12 and string-inspired theories [13] are studied extensively. Recently a 
theory of /(T) gravity has been received attention. Models based on modified teleparallel gravity were 
presented, in one hand, as an alternative to inflationary models |141 I15j . and on the other hand, as an 
alternative to dark energy models |16j . New spherically symmetric solutions of black holes and wormholes 
are obtained with a constant torsion and for the cases for which the radial pressure is proportional to 
a real constant, to some algebraic functions f(T) and their derivatives, or vanishing identically [21| . In 
the same way, an algebraic function /(T) is obtained through the reconstruction method for two cases 
and the study of a polytropic model for the stellar structure is developed [22]. Moreover, /(T) gravity is 
reconstructed according to holographic dark energy is explicitly presented in [23^ and latter an anisotropic 
fluid for a set of non-diagonal tetrads in /(T) gravity explored generating various classes of new black hole 
and wormhole solutions 24 . Also, many works have been done in order to check whether /(T) gravity 
can present results consistent with the many advances in cosmology and astrophysics 25[. Recently, 
Bamba et al investigate the reconstruction of power law model, exponential model and logarithmic model, 
able to reproduce some of the future flnite time singularities and also discuss the thermodynamics near 
these singularities 261 (For other works about future flnite time singularities, see [27]). Also, future 
singularities vifith the presence of a viscous fluid as well as other interesting properties, have 
been already studied (in the context of General Relativity and f(R) gravity) (28] . 

In the present paper we investigate the /(T) gravity models that are able to reproduce the four 
types of flnite time future singularities from a suitable choice of Hubble parameter. A parameter /3 in 
the expression of the Hubble parameter plays an important role in specifying these singularities. The 
algebraic function f{T) is assumed as the sum of the teleparallel term T and an algebraic function g{T) 
with which all the task is done. According to some values of the parameter /3, differential equation of 
g{T) are established and solved in some ways. The algebraic function f{T) for each type of singularity 
is obtained from each expression of g{T). 

On the other hand, we notify that the presence of finite-time future singularities may cause serious 
problems in the black holes or stellar astrophysics [29] . A way to probe the possible avoidance of these 
singularities is considering that the fluid possesses viscosity. This is the second purpose of this work. As 



previously mentioned, the models that lead to future finite time singularities have been reconstructed 
considering a non viscous fluid. Thus, we introduce the viscosity and investigate its effects on the 
singularities. We see that when the constant viscosity or the viscosity proportional to \J —T is considered, 
in general, the singularities are robust against the viscosity. However, when the viscosity is proportional 
to (— T)"/^, for some values of the parameter n, the viscosity may cure the Big Rip and the Big Freeze. 
The paper is organized as follows. In Sec. 2, the /(T) gravity formalism and the field equations 
are presented. The Sec. 3 is devoted to the presentation of the Hubble parameter, the classification of 
future finite time singularities. Suitable scale factor coming from the Hubble parameter are presented 
according to the values of the parameter /3 for obtaining the algebraic /(T) function. The viscosity is 
introduced in the Sec. 4 and its effect is investigated as the singularities are approached. The conclusion 
and perspectives are presented in Sec. 5. 

2 /(r) gravity and field equations 

Let us define the notation of the Latin subscript as those related to the tetrad fields, and the Greek 
one related to the spacetime coordinates. For a general specetime metric, we can define the line element 
as 

ds^ = g^ydx^dx'' . (1) 

The projection of this line element can be described in the tangent space to the spacetime through the 
matrix called tetrad as follows: 

ds^ = g^^dxf'dx" = riijO'e^ , (2) 

dx^^Ci^e' ,e'^e\dx'' , (3) 

where r]ij is the metric on Minkowski's spacetime and e^^e^ ^ = (5^, or e^^e^ = 51- 

The action for the theory of modified gravity based on a modification of the teleparallel equivalent of 
General Relativity, namely f{T) theory of gravity, coupled with matter Lm is given by [15 1 [TS t fTO l [20 ] 



S^^^jd\e[T + g{T) + L„,]^ [A) 

where e — det{e^ ) = \J~-g. Here, G is the gravitational constant and c the speed of the light. From now, 
we will use the units SttG = c = 1 . The teleparallel Lagrangian T is defined as follows 

where 

T\.=e,\d^e\-d,e\\ (6) 

S, ^^ = \{K^% + S^pT'\ - 5;t\), (7) 



and Kf^"^ is the contorsion tensor 

K^% = -\{T^''p-T'^^p-T^n- (8) 

The field equations are obtained by varying the action with respect to vierbein ej^ as fohows 

- e-'d,ieS, ''"Kl + gr) - e,'T^^,S^ ''^gr + S, ^'"d.^grT - \eW + g{T)) = ef% ^ (9) 

where gx = 9'{T) and gxT = g"{T) and T the energy momentum tensor. Now, we take the usual 
spatially-flat metric of Friedmann-Robertson- Walker (FRW) universe, in agreement with observations 

3 

ds^ = dt^ - a(t)2 ^(da;*)^ (10) 

where a{t) is the scale factor as a one-parameter function of the cosmological time t. 

Let us assume first that the background is a non- viscous fluid. Using the Friedmann-Robertson- Walker 
metric and the perfect fluid matter in the Lagrangian ^ and the field equations ^, one obtains 

T = -6iy2 , (11) 

m^ = Peff , (12) 

-3i7' -2H = p,ff , (13) 

where Pe// and Peff denote respectively the effective energy density and pressure of the universe and 
defined by 

Peff = P-^a- 6H^9T , (14) 

Peff - p+^g + 2(3H^ + H^gT-2AHH^gTT , (15) 

where H is the Hubble parameter and defined hy H — a/ a. Using ([T4|) and (fTS]) . and combining ([72]) and 
(fT5)l . one gets 

2T5rr + 5T + ^^^4-^ + 1 -0 , (16) 

2H 

where we used the barotropic equation of state p — up. Then, for a given scale factor corresponding 

to a future finite time singularity, the action may explicitly be reconstructed by solving the differential 

equation (fT6|) . 



3 Future finite time singularities 

We propose to find in f{T) gravity, models that reproduce the four types of finite time future singu- 
larities from the Hubble parameter [35] 

H^h{ts-ty^ + C , (17) 



where h, C and t^ are positive constants and t < tg. These constraints are imposed to the parameter for 
providing an expanding universe. The parameter /? can be a positive constant or a negative non-integer 
number. Then, as the singularity time ts is approached, H or some of its derivatives and therefore, the 
torsion, diverge. C is essentially relevant near the singularity only when /3 < (where we denote it as 
C = Hg, the Hubble parameter at the singularity time), and then, we can assume C = Q when /3 > 0. 
The finite-time singularities are classified in the following way [301 IS] 

• Type I (Big Rip): for i — > t^, a — > oo, pe// -^ co and \peff\ -^ co at t — tg. This corresponds to 
13 = I a.nd /3 > 1. 

• Type 11 (Sudden): for f — > t^, a — > as, Peff — ^ Ps and |pe//| — > oo. It corresponds to —1 < /? < 0. 

• Type III (Big Freeze): for i ^- ts, a —> as, peff — > oo and \peff\ — > oo. This corresponds to 
0<^<1. 

• Type IV (Big Brake): for t ^ tg, a ^ gq, Peff -^ 0, Pejf —J- and higher derivatives of H diverge, 
which corresponds to the case /3 < — 1 but /3 is not any integer number. 

Let us now investigate the /(T) gravity models for which the finite time future singularities could 
occur, when ([T7| is assumed. 

3.1 Big Rip singularity models without viscosity 

This sort of singularity may appear for /? = 1 and /3 > 1. Let us treat the cases separately. 

The case /3 = 1 

In this case, the corresponding scale factor can be written as 

a{t) = ao {ts - ty^ , (18) 

and then, the first derivative of the Hubble parameter reads 



H = h{ts~t) ^ . (19) 



Hence, Eq. ([T6| becomes 



po{i+u.){ts~tr^'^-^+' 



2TgrT + 9T + ^ ^^^^^ + 1 = 0. (20) 



In the other hand, using pT|). with the scale factor p^ . one can write 



'^-'-(-i^) ''' • (21) 



Thus, Eq. (PU)) takes a new form as 






The general solution of ([22l) reads 
-1 



9iT) 



1 + 2B 



1 + 2B 



A 



T 



l + B \ 6/i2 



2(l + 2B)Ci\/^ 



T + C2 



where Ci and C2 are integration constants, and A and B defined respectively as 



A = 



Po(l+^) 

3(l+cj) 




2hd, 



B = — h{l+uj)-l 



The corresponding algebraic function /(T) reads 



!{T) 



AT 



T 



-— -2CiT^/-T + C2 



(23) 



(24) 



(25) 



{l + 2B){l + B) V 6/i2^ 

Initial condition may be applied for finding the respective values of the constants Ci and C2 ■ We can 
follow the same process as in [M]. We assume that at the early time, that we denote tg, the corresponding 
value (the initial value) of the torsion scalar is Tq such that 



dT 



= -12h^ - 



t = to 






3/2 



The initial conditions imposed to the functions / read 



if)t=tn - ^0 , 



dT 



t=to 



Making use of these initial conditions, ([?f|) . one gets 

B' 



Ci 



1 



sv^n 



1 



A 



1 + 2B \ 6h^J 



, C2- 



2Tn 



1 + 



A{1 - 2B) 



Tn 



2{l + B){l + 2B) V 6/i2 



(26) 



(27) 



(28) 



Then, the algebraic function (1^51) . with the constants (P5|) . leads to the Big Rip when the fiuid is free of 

viscosity. 



The case (3 > 1 

In this case, the corresponding expression of the scale factor is 

a{t) = age 3^^i , 
and the first derivative of the Hubble parameter reads 



H ^ Ph (ts - t) 



-/3-1 



Thus, the Eq. p^ becomes 



2TgTT + 5t 



^-^ ■,,. . ^ e 1^^^ + 1 



2phal^'+'^'> 



= 



(29) 



(30) 



(31) 



Using Eq. (fTTI) and the scale factor (f29| . one gets 






and the differential equation pip can be written as 



Note that this equation is more complicated and cannot be solved trivially. However, it can be solved as 
the singularity is approached, i.e., when t — >■ i^. Then, as the singularity is approached, the trace diverges 
(T — > —00), and ([33| becomes 

2TgTT + 9T + 1 = , (34) 

whose solution is 

5(T) = -r + 2C3\/^+C4 , (35) 

where C3 and C4 are integration constants. Consequently, f(T) is written as 

/(T) = 2C3\/^+C4 . (36) 

In this case, the initial value of the first derivative of the torsion with respect to the cosmic time reads 

Making use of ([57]) and the initial conditions, (P?]). one gets 



C3 = -V^ , C4 - -To . (38) 

Then the corresponding algebraic f{T) is 



fiT) ^-2y/nf -To . (39) 

This model also leads to the Big Rip in a universe dominated by a non-viscous fluid. 

3.2 The Big Freeze models without viscosity 

The Big Freeze appears for < /? < 1, and the corresponding scale factor is the same as in ([2^ . 
which leads to (|33)) . However, as we are dealing with the singularity of type IH, Eq. (|33|) becomes 



whose general solution is 

9{T) = -T 






h{/3 - l)a^' 
where C5 and Cg are integration constants. Then, /(T) is written as 



(41) 



u^a T\ 3(l+w) " 



6/1= 



(42) 



By using I^E^ and the initial conditions (P7|. the constant C5 and Cg are determined as 



a 



Po(l+^)V6 f_T^ 
2 3(1+") I 6/^2 



, Cq — Tq 



oo{l + uj) (3+1 



^/,(^_l)„3(i+-) /3 



Then, the algebraic function (j42|) is characteristic of the singularity of type III. 






(43) 



3.3 The Sudden singularity models without viscosity 

In this case, — 1 < /3 < 0, and the corresponding scale factor reads 



-(t.-o(g.- "vr' 



a(i)=aoe V ^-^ j ^ (44) 

and the first derivative of the Hubble parameter remains the same as in ([50]) , and Eq. p^ becomes 



po(l + u;)(t,-t)^+i 3(i+^)(t.-t)(g.^ '"%-jr" ) 



By using ([T7)) . with C substituted by Hg, as explained above and pT|) . one obtains 



(45) 



(i. - = 



6/i2 /i 



(46) 



where T, denotes the torsion scalar at the singularity time, and may not be infinity. Using (|46p . Ea. (j45p 
becomes 



2TgTT + gr 



Po(l+^) 



6/i2 h 






,3(l+-)[/^-^]"^(ff-^[/^-^]) +1 = 







(47) 



This equation is not trivial and cannot be solved analytically. However, as the singularity time is ap- 
proached, with — 1 < /3 < 0, it becomes 



2TgTT + gr 



Po(l+^) 



T Hs 
6/i2 h 



+ 1 = , 



(48) 



whose general solution reads 



9{T) 



2-^lihA 
/3-1 



-4T Hs 

2 - 

9/i2 h 



y/^6T-6H, 



2C7V^+Cs 



(49) 



where C7 and Cg are integration constants, and A previously defined in (|24[) . The algebraic function 
/(T) is then written as 



fiT) = 



3+1 
2^^/3hA 

/3-1 



2 ^ 

9/i2 h 



V^6T-6Hs 



+ 2C7V^ + Cg 



(50) 



The constants C7 and Cg can be determined in the same way as in the previously cases. In this case, the 
first derivative of the torsion scalar with respect to the time, at initial time, is written as 



'^Lr-M-'i 



and (P7| is valid. Then, we obtain the constants 



1/2 



To Hs 



C7 = - 



5/9 + 2 

2^i^l3A 
V3(/3 ~ 1) 



V^ + 



I3A{2 + V6)2 



3 + 1 



Cs — To 



3+1 
2T9-/3/iA 

/3-1 



-4To g, 
9/i2 /i 



v/=6li;-6i/, -2C7v/=7i^ 



(51) 

(52) 
(53) 



Then, the algebraic function ([SO]) . with the constants C7 and Cg respectively in (15^ and (15^ . is the /(T) 
models that may produce the singularity of type II. 



3.4 The Big Brake models without viscosity 

Here, the scale factor is the same as in the case of the Sudden singularity, just that in this case, 
/? < — 1. Then, it can be easily observed that the differential equation (|48p is also valid in this case. 
Consequently, the algebraic function (|50|) can lead to the Big Brake. The difference which appears here, 
with respect to the sudden singularity, is the values of the parameter /3. 



4 Analysing the possible avoidance of the singularities in a vis- 
cous fluid 



Let us assume the fluid equation of state in the following form 



(54) 



where C{p) is the bulk viscosity and in general depends on p. According to the thermodynamics grounds, 
the quantity ({p) has to be positive in order to guarantee the positive sign of the entropy change in an 



irreversible process. In this case, the stress-energy tensor T^^ is written as 

%i, = pu^u„ - [up - 3HC{p)] {gp^i, - u^u,,) . (55) 

The equations of motion ([T^ and p^ can be written as 

3H^ = p + PT , (56) 

-2H-3H'^ = p + PT , (57) 

where the modified gravity part is formaUy inchided into the modified energy density pr and the modified 
pressure pr as foUows 

PT = -^g-6i/V , (58) 

Pt = ^g + 2(3H^ + H)gT-24HH^gTT ■ (59) 

If we assume that the ordinary and the dark fluids of the universe do not interact, the equation of 
continuity related to the ordinary viscous fluid reads 

p + 3Hp{l+w)=9H^C{p) ■ (60) 

As we have seen in the section 3, when the fluid does not possess viscosity, the four type of future 
finite-time singularities may appear. Now by introducing the viscosity, we analyse its effect near the 
singularities. Let us start with the constant viscosity case. 

4.1 The constant viscosity case: ({p) = (q 

In order to analyse the effect of the viscosity on the Big Rip, Sudden, Big Freeze and Big Brake singular 
models, it is worth considering the behaviour of the conservation law in Eq. (1601) near the singularities. 
By using Eq. p7)) . one can write Eq. (pO)) as 

p + 3p (I + w) [its - ty^ + C] ~ 9Co/i' (ts - t)~^^ + ISCoC/i (ts - t)-^^ + XqC . (61) 

Solutions of Eq. (ICTI) can be found according to different values of the parameter (3 as 

9h\o 



Pvis 
Pvis 
Pvis 
Pvis 



{ts-t){l + 3h + 3huj) ' 
3feCo 

(2/3-l)(i,-t)2/3-i ' 

9hHsCo SHsCo 

(/3-l)(i,-i)/3-i + 1 + c^ 
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for /3 = 1 , 


(62) 


for /3 > 1 , 


(63) 


for < /3 < 1 , 


(64) 


r /3 < , iJ, 7^ , 


(65) 



where we used C = for the Big Rip and the Big Freeze as previously discussed, since this constant is 
not relevant in these cases, while for the Sudden and the Big Brake, the constant C is taken equal to the 
Hubble parameter Hg at singularity time. Since Hs plays a crucial role in these two cases, it has to be 
different from zero. 



• 



For /? = 1, (Big rip), the energy density of the viscous fluid behaves as (tg — t)^ , while the Hubble 
parameter behaves as {tg — i)^^. We see that the viscosity part of the energy density diverges more slowly 
than the Hubble parameter. Hence, in this case, the constant viscosity cannot avoid the occurrence of 
the Big Rip for /? = 1. 

• For /3 > 1, the Hubble parameter behaves as (t^ — t)^"^^ , while the energy density from the viscous 
fluid behaves as {ts —t)^^, then, is less that the Hubble parameter. Hence, also in this case, the Big Rip 
cannot be avoided by a fluid with constant viscosity. 

• For < /3 < I, the energy density from viscous fluid behaves as {tg — i)^^^'^, while the Hubble 
parameter behaves as (ts — i)~^'^. We see in this case that the energy density of the viscous fluid diverges 
less that the Hubble parameter. Hence, the Big Freeze cannot be avoided by constant viscosity. 

• For — 1 < /3 < 0, both the Hubble parameter and the energy density from the viscous fluid are 
finite. The task here is to consider the behaviour of the first derivative of the Hubble parameter H 
and the pressure of the viscous fiuid. Note that as the Sudden singularity is approached, H behaves 
as {ts — t)^^^^ , while the pressure Pms of the viscous fluid behaves as {ts — t)~^^ . It appears that the 
pressure is finite while H diverges. Then, the Sudden singularity is robust against the constant viscosity. 

• Since H diverges in this case of constant viscosity, the higher derivatives of H also diverge, then the 
Big Brake cannot be avoided by constant viscosity. 

We conclude that in /(T) gravity and with the equation of state ([54| for a fiuid with constant viscosity, 
all the singularities are robust against the viscosity. 



4.2 The viscosity proportional to \/—T 

From Eq. pT|) one can observe that this case is that in which the viscosity is proportional to the 



Hubble parameter. Let us then consider C = ta/— 6T/2 ~ 'StH, where r is a positive constant. The 
different solutions of (|60|) according to the values of the parameter /3 read 

'- = {ts~tn2 + lh + 3hco) ' '°^ ^ = '' ^^^^ 

P- " {l+.){ts-tfP ^ ^°'' ^>^' (^^) 

27h!^T , , 

P- " (3/3-l)(t.-t)3^-i ' for 0</3<l, (68) 
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Let us now discuss the effects of the viscosity near each type of singularity. 

• For 13 = 1, the energy density of the viscous fluid diverges like H^. For small values of r, one 
can conclude that the Big Rip cannot be avoided. However, for large values of r the Big Rip could be 
prevented. 

• For /3 > 1, pyis diverges like H^. Thus, for small values of r, the Big Rip cannot be avoided, while 
for large values of r it can be avoided. 

• For < ^ < 1, py,s diverges like {t^ - tf-'^l^ while the H'^ diverges like (i^ - 1)-"^!^ . For < ;3 < 1/3, 
Pyis is finite while for l/3</3<l, it diverges. However, in the case pms diverges, it diverges less than 
H^ . Thus the Big Freeze is robust against the viscosity. 

• For — 1 < /3 < 0, pvis tends to 9tHs/{1 + w), then, the pressure Pms tends to 9(jJtHs/{1 + w). At 
the same time the first derivative of the Hubble parameter, ij, diverges. Thus, for the small values of 
r, the viscous fiuid cannot influence the Sudden singularity. However, for large values of r, the viscous 
pressure could dominate over the H end then, the Sudden singularity may be avoideqfl. 

• For /3 < — 1, one can analyse the behaviour of the second derivative of the Hubble parameter, H , 

1/2 

near the Big brake and compare it with the behaviour of the second derivative of pj^^ . One can observe 

1/2 

that H diverges only if —2 < /3 < — 1. For the same interval, pj^^ also diverges as the singularity is 
approached, but less than H . Hence, the viscous fluid in this case cannot influence the Big Brake. 

4.3 More general case: C, proportional to (— T)"/^ 

In this general case, we consider the bulk viscosity as 

C = 'r(^^y =r{3Hr , (70) 

where n is a natural number, and r a non null positive constant. Hence, the energy conservation leads to 

p + 3ff(a; + l) = 9H2(3iJ)"r , (71) 

from which we obtain the behaviour of the energy density of the viscous fluid as 

'^^ = ^i^oin^ A ts-tr-' , for (3 = 1, (72) 

^^^)\^t^_ty(n+i)0 ^ for /3>1 , (73) 



Pvi 
Pvi 



a; + 1 
(3/i)"+2 



(2 + 71)^-1 



— (i, - t)i-^("+2) , for < /? < 1 , (74) 



CiHsT+^hT. ,i_. , (3iJ,)"+i 



T 



Hs[p - 1) l+w 

• For (3 = 1, H^ diverges like {ts — t)^"^ , while pyis diverges like (t^ — t)^^^". Here the situation 
becomes discussable due to the presence on the parameter n. Hence, one observes that for < n < 1, as 



^This is possible only if r is sufficiently large for dominate over the divergence of H. 
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the singularity is approached, H^ diverges more than p-uis- Then, the Big Rip cannot be avoided by the 
viscous fluid. However, for n > 1, the pyis diverges more than H^ and then, the Big Rip can be avoided. 

• When /? > 1, Pvis diverges like {tg — t)^^"+^)^, while H^ diverges like (i^ — t)^^- Here, we see that 
for n < — 1 + 2/(3 (with 1 < (3 < 2), H^ diverges more than p^s. Hence, the Big Rip is robust against 
the viscous fluid. But when n > 2/3—1, p„is dominates over H^. Hence, the Big Rip may be avoided 
from the effect of the viscosity. 

• For < /3 < 1, iJ2 diverges like {ts-ty^, while p^^ diverges like (is-i)^~'^'"+^^. When n < -2+3/(3 
(with 1 < /3 < 3), the energy density of the viscous fluid diverges less than H^ and the Big Freeze cannot 
be avoided. But for n > —2 + 3//3, the energy density of the viscous fluid dominates over the background 
and the avoidance of the Big Freeze becomes possible. 

• For — 1 < /3 < 0, both the viscous energy density and H^ are finites. In this case, it is necessary 
to compare the behaviour of the pressure of the viscous fluid with the first derivative of the Hubble 
parameter, i.e. H. Note that the pressure of the viscous fiuid is also finite but depends strongly on the 
parameter n, that is, Pms = (^t{3Hs)"^^ /{i^ + 1), while H diverges like {ts — t)^^'^. Then, a priori, the 
Sudden singularity cannot be avoided by the viscous fluid. However for large values of n, Hs and r, the 
pressure of the viscous fluid may dominate over H and the Sudden singularity could be avoided. 

• For /3 < — 1, if the values of the parameters n, Hs and r are very large, the viscous fluid can influence 
the feature of the Big Brake, but not necessary avoid they. 

As conclusion, we see that for a viscous fluid those viscosity is proportional to (— T)"/^, with the 
equation of state ([Ml) , the Big Rip and the Big Freeze can be eliminated for some values of n. However, 
in the case of the Sudden and the Big Brake, for large values of n, Hg and r, the viscous fluid influences 
the singularities but does not necessary avoid they. 

5 Conclusion 

We considered in this work the modified teleparallel theory, known as f{T) theory, where T is the 
torsion scalar. In a specific way, the algebraic function /(T) is taken as the teleparallel term T plus 
the algebraic function g{T). The equation of motion of the theory is used and differential equation is 
established with the algebraic g{T). The expression ()17|) is assumed for the Hubble parameter where the 
parameter (3 plays an important role in specifying the type of singularity. Besides the parameter (3 which 
is sufficient for characterizing the singularities, the Big Rip and the Big Freeze, we need to introduce 
the parameter Hs, which substitutes the constant C, for specifying the Sudden singularity and the Big 
Brake. Then, the differential equation is solved in each case and the corresponding algebraic function 
/(T) which may lead to each type of singularity is obtained. All this is done considering a fluid without 
the viscosity. 
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In order to probe the possible avoidance of the singularities, we introduce the bulk viscosity (^(p) 
in three ways. The first case is when the viscosity is constant and then we observe that in general, the 
viscosity is inefficient against to the singularities. In the second second where the viscosity is proportional 
to \/— T, we see that for small values of the parameter r the singularities are robust again the viscosity, 
while for large values of the this parameter, just the Big Rip and the Big Freeze may be cured by the 
viscosity (the Sudden and the Big Brake remain robust in this case) . In the third case where the viscosity 
is proportional to (— T)"/^, we observe that independently of the values of r, the Big Rip and the Big 
Freeze may be avoided for some values of the parameter n. However, the Sudden and the Big Brake could 
be cured only for large values of n, r and Hg. 
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